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We have systematically studied the mechanical stability of all noble metal carbides with the rock- 
salt structure by calculating their elastic constants within the density function theory scheme. It 
was found that only four carbides (RuC, PdC, AgC and PtC) are mechanically stable. In particular, 
we have shown that RuC, PdC, and PtC have very high bulk modulus, which has been remarkably 
observed by the most recent experiment for the case of PtC. From the calculated density of states, 
we can conclude that these compounds are metallic, like the conventional group IV and group V 
transition metal carbides. 
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Transition metal carbides (TMCs) have attracted 
much attention due to their excellent physical, chemical 
and mechanical properties. fl In the past, most work 
was focused on the group IV and group V TMCs,0, 0, 
S while noble metal carbides have been only stud- 
ied spectroscopically by few researchers |j[ I9I Hfjj because 
the crystalline samples are hard to obtain. Very recently, 
shortly after the synthesis of a novel noble-metal nitride, 
PtN,[ll| Shigeaki Ono et al. 12] have succeeded for the 
first time in synthesizing platinum carbide (PtC) at high- 
pressure and high-temperature by the laser-heated dia- 
mond anvil cell technique. According to their report, the 
new PtC has a rock-salt (RS) structure, with a high bulk 
modulus of 301 (±15) GPa. 

While the structual and electronic properties of earlier- 
synthesized PtN have been extensively studied in recent 
theoretical[llJp EEEESElEllSlIllilliiand 

experimental|2J] work, very few attention has been payed 
to the PtC system. Recently, Li et aZ.jl^] have given a 
first-principles calculation of PtC. They preferred PtC 
crytallized in zinc-blende (ZB) structure since their re- 
sults showed that the RS-PtC is mechanically unstable. 
In contrast, after a systematic calculation, and based on 
two separate DFT schems, we found that the RS-PtC is 
mechanically stable 26] and the compressibility behavior 
of RS-PtC is more comparable with that of the exper- 
iment. To illustrate, here we show in Fig. 1 the cal- 
culated pressure-volume relation for both ZB- and RS- 
PtC. For comparison, the experimental data^l are also 
plotted. Clearly, one can see that the calculated com- 
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FIG. 1: Calculated Pressure- Volume data (average LDA and 
GGA results) of PtC with RS and ZB structure, PtC 2 with 
Pyrite structure and experimental data. Dashed curves are 
the third-order Birch-Murnaghan equation fit with B=289 
(239, 257, 301) GPa and B'=5.Q (4.7, 4.6 5.3) for the RS 
(ZB, Pyrite, exp.). 

pressibility ratio for RS-structure is in good agreement 
with the experimental result while it is not so for ZB- 
structure. Therefore, our results firmly show that the 
observed phase of PtC under high- pressure is rocksalt, 
although it is difficult to distinguish the structures of 
RS from ZB only by X-ray diffraction measurement |27j . 
We have also found 26] that the total energy of the ZB 
structure is lower than that of the RS structure (about 
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1 eV/unit cell). From this sense, the synthesized PtC 
in Ref.[12] is metastable. In addition, we have further 
calculated the elastic constants of PtC2 with the flu- 
orite and pyrite structures, since the former structure 
was once proposed to explain the high bulk modulus of 
platinum nitride and the latter structure has been re- 
cently observed |24j in platinum nitride. Our calculation 
revealed that the fluorite PtC2 is mechanically unsta- 
ble. On the other hand, the pyrite PtC2 is calculated 
to be mechanically stable but with a low bulk modu- 
lus (B=236 GPa, B'=4.7), much lower than that of the 
experiment, [l^ What's more, there also exsits a big dis- 
crepancy between the calculated pressure- volume data of 
pyrite PtC2 and those of experiments (see Fig.l). There- 
fore, we excluded the possibility of existence of zinc- 
blende or pyrite structure in the experiment. |l2| 

To appreciate the consistency of our calculation with 
the experimental observation on PtC and its very high 
bulk modulus, in this paper we have gone beyond to 
search in a systematic way for other possible mechani- 
cally stable noble-metal carbides with the RS structure 
by calculating the relevant elastic constants. Surpris- 
ingly, we have found that RuC has even higher bulk 
modulus of ^370 GPa. Among the whole noble metal 
(Ru, Rh, Pd, Ag, Os, Ir, Pt, and Au), Ag and Pd were 
also found to form mechanically stable carbides. To- 
gether with RuC and PtC, they form a new family of low 
compressible compounds. The partial density of states 
(PDOS) of RuC, PdC, AgC, and PtC were calculated, 
which revealed that the hybridization of metal 5c?-like 
and C 2p-like states is responsible for the strong covalent 
bonding of these compounds. The PDOS also revealed 
that these compounds are of metallic character. 

The first-principles calculations were performed by the 
plane- wave pseudo-potential (PW-PP) method as imple- 
mented by CASTEP code[28|. The ultrasoft pseudopo- 
tential (USPP)0| was used to describe the interaction 
between ions and electrons. The exchange and correla- 
tion terms were described by the local-density approxi- 
mation (LD A) [ 30| and generalized gradient approxima- 
tion (GGA)[3l|. For the Brillouin-zone sampling, the 
Monkhorst-Pack (MP) scheme^ with a mesh of (12, 
12,12) was adopted. We chose Ecut =550 eV for all the 
calculations. It was tested that with even more strict 
parameters the total energy can be converged within 
0.003 eV/atom for all the cared systems. For the me- 
chanical stable phases, the elastic constants were recalcu- 
lated with the full-potential linearized augmented plane- 
wave (LAPW) method^ (implemented by WIEN2K), 
wherein a non-overlapping muffin-tin (MT) sphere radius 
of 2.2 (2.1, 2.3, 2.3) a.u. for Ru (Pd, Ag, Pt) and 1.7 (1.5, 
1.7, 1.7) a.u. for C atoms in RuC (PdC, AgC, PtC) com- 
pound was employed. The valence wave functions inside 
the MT spheres were expanded into spherical harmonics 
up to Z=10 and the potential up to l=A for the Ad carbides 
and 1=6 for the 5c? carbide (PtC). The energy that sepa- 



rates the core and valence states was set to be -6.0 R y for 
the Ad carbides and -7.0 R y for the 5d carbides. Here, an 
APW+lo[H type basis set for certain electronic states 
was also used. We used 5000 k points for sampling the 
Brillouin zone and R m t^<K max was taken as 8.0 (where 
R m t is the muffin-tin radius and K max is the plane-wave 
cutoff). A fully relativistic calculation was performed for 
core states, whereas the valence states were treated in a 
scalar relativistic scheme. 

After getting the equilibrium geometry configures, in 
the PW-PP method we applied the so-called "stress- 
strain" method to obtain the elastic constants. For the 
cubic crystal, there are only three different symmetry 
elements (en, C12 and C44). A single strain with non- 
zero first and fourth components (xx and yz) can give 
stresses relating to all three of these coefficients, yielding 
a very efficient method for obtaining elastic constants for 
the cubic system. The elastic constants can be obtained 
from the stress-strain relations and then the bulk mod- 
ulus is obtained from the elastic constants by equation 
5=(cii+2ci2)/3. In the LAPW method another scheme 
was adopted (For details read document in WIEN2K 
package). In this scheme, three types of deformation were 
performed to produce three linear equations for en, C12, 
and C44. 

Table I lists the calculated clastic properties of the me- 
chanically stable noble metal carbides. One can see that 
the results from the PW-PP method and those from the 
LAPW method fit very well. An exception is the ampli- 
tude of C44 in PtC and RuC; the inconsistency between 
the two methods remains unclear. However, since the 
amplitude of C44 does not affect the value of the bulk 
modulus of these compounds, we do not care about it at 
present. As it is well known that in first- principles calcu- 
lations the LDA (GGA) usually underestimates (overes- 
timates) the lattice constant and overestimates (under- 
estimates) the bulk modulus, hence we adopt in the PW- 
PP method to use the average of LDA and GGA values 
as the re-scaled estimates. The responses of the crystal 
under hydrostatic pressure, rhombohedral, and tetrago- 
nal distortions are measured by the bulk modulus B, the 
normal shear modulus G' = C44, and the tetragonal shear 
modulus G" ={c\\-c\2) jl, respectively. The mechanical 
stability of a crystal requires the strain energy to be pos- 
itive, which for cubic crystal implies pjif 

C44 > 0,Cn > |ci2 1, en + 2ci2 > 0. (1) 

From Table I we can see that the listed four carbides are 
all stable because their elastic constants satisfy formula 
Q). One can find that among these compounds RuC 
has the highest bulk modulus with the value ~347 GPa, 
which is comparable with that of cubic BN (367 GPa 36] ). 
The PtC has a theoretical bulk modulus of 288.7 GPa, 
in good agreement with the experimental value of 301 
GPa. Thus Table I shows that the rocksalt RuC, PdC, 
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TABLE I: The calculated equilibrium lattice parameters a(A), elastic constants (GPa), bulk modulus B (GPa), tetragonal 
shear modulus G" (GPa), poly crystalline shear modulus G=(cii-d2+3c44)/5 (GPa), Young's modulus E (GPa) and Poisson's 
ratio of four mechanically stable noble metal carbides. 
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together with PtC, form a new family of low compressible 
compounds. 

We notice that the trends of bulk modulus in these 
carbides are in some way related to the number of va- 
lence electrons. This can be seen by the fact that the 
value of bulk modulus decreases with increasing the to- 
tal number of valence electrons (12 per primitive cell for 
RuC, 14 for PdC and PtC, and 15 for AgC). To further 
investigate the origin of the bulk modulus, we have cal- 
culated the partial DOS of these carbides. The results 
are shown in Figs. 2. For RuC (Fig. 2a), there ex- 
ists a minimum in DOS near the Fermi level, giving the 
value of about 0.6 states/eV per primitive cell. The states 
between -15 eV and -11 eV are mainly composed of C 
(2s) states with a little contribution from Ru (Ad) states. 
There exist two peaks of the Ru (Ad) states between -9 
eV and the Fermi level. The strong hybridized states of 
Ru (4d) and C (2p) states between -9 eV and -A eV is 
responsible for the covalent bonding of RuC. For PdC 
and PtC which contains two more valence electrons than 
RuC, the extra electron fills the anti-bonding d orbital 
of transition metal atoms, thereby decreasing the bond- 
ing strength of such compounds. As a consequence, the 
bulk modulus of PdC and PtC are much lower than that 
of RuC. The same rule can apply to AgC. One may ask 
the question why there is an invisible difference of bulk 
modulus between PdC and PtC. This can be explained 
as follows: the presence of 4/ electrons in the latter com- 
pound, which should repel (i-orbitals out of core regions 
thus enhance the d-p bonding (see the DOS patterns of 
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FIG. 2: Partial DOS of (a) RuC, (b) PdC, (c) AgC and (d) 
PtC with rock-salt structure. 



Fig. 2b and Fig. 2c), giving a higher value of bulk mod- 
ulus for PtC. In another way for a quantitative analysis, 
we have investigated the Fermi energy position of PtC, 
PdC and AgC compared to RuC, whose Fermi energy lo- 
cates at d states of transition metal atoms. We find the 
more the Fermi level separated from d states, the lower 
the bulk modulus approaches. The trend is compatible 
with above proposed mechanism for bulk modulus de- 
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creasement. Furthermore, from Figs. 2 one can conclude 
that all these four carbides are metallic. 

As mentioned above, the tetragonal shear modulus 
G" measures the response of a crystal to the volume- 
conserving tetragonal shear strain. It is thought that the 
stretching of metal-C bonds and the bending of metal- 
metal bonds are involved in such procedure |37j. Because 
en is determined by the nearest-neighbor interaction as 
well as bulk modulus B, thus when cn»c\2, G" have 
the same trend as B, as we observed in RuC and PdC 
in Table I. However, when en is comparable to C\2, the 
shear modulus trend does not need to be similar as that 
of bulk modulus. For example, in the present calcula- 
tions, the bulk modulus of PtC is 288.7, much larger 
than that of AgC (158.8 GPa). However, its G" (27.8 
GPa) is smaller than that of AgC (42.7 GPa). 

The nature of hardness has been extensively investi- 
gated and many new models have been proposedJH, US 
14(1 l4ll |42| . Although for the strong covalent materials, 
the hardness can be directly evaluatedj^ 0] > f° r the 
partially covalent transition-metal based hard materials, 
the polycrystalline shear modulus G [(en — ci2 + 3c44)/5] 
is still considered as the most important parameter be- 
sides the bulk modulus, governing the indentation hard- 
ness. Our calculated polycrystalline shear modulus G 
of RuC is 65.2 GPa, which is the highest value of these 
four noble metal carbides. This value is just comparable 
to that of pure Pt (61 GPa0), remarkably contrasting 
the very high value of bulk modulus. It was also found 
that the highest shear modulus C44 of these four com- 
pounds (PtC: 49.4 GPa) is much lower than that of con- 
ventional materials. For example, the C44 of pure Pt is 
76.4 GPa. 0] Therefore these four potential noble metal 
carbides are not expected to withstand shear strain to a 
large extend and may have very low hardness. 

In conclusion, we have performed systematic ab initio 
total energy calculations on the elastic stability of all the 
rock-salt structured noble metal carbides. It was found 
that four carbides, namely RuC, PdC, AgC and PtC, 
are mechanically stable. Most of these potential transi- 
tion metal carbides have low compressibility though their 
hardness is not expected to be high. The low compress- 
ibility was found to be related to the occupancy of the 
valence states. From the calculated electronic structure, 
all these carbides is metallic. Since the rock-slat PtC has 
been recently synthesized under high pressure and high 
temperature, we expect the other three carbides can also 
be prepared under certain extreme conditions. 
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